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Optical lattice loaded with cold atoms can exhibit a tunable photonic band gap for a weak probe 
field under the conditions of electromagnetically induced transparency. This system possesses a 
number of advantageous properties, including reduced relaxation of Raman coherence and the asso- 
ciated probe absorption, and simultaneous enhancement of the index modulation and the resulting 
reflectivity of the medium. This flexible system has a potential to serve as a testbed of various 
designs for the linear and nonlinear photonic band gap materials at a very low light level and can 
be employed for realizing deterministic entanglement between weak quantum fields 
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I. INTRODUCTION 

The properties of waves in spatially periodic media are 
studied in various branches of physics, including acous- 
tics, electromagnetism and quantum mechanics, initially 
called wave mechanics. Some of the fundamental prop- 
erties of solid-state systems stem from the interplay of 
wave-like behavior of amplitudes describing the mobile 
electrons and the spatially periodic potential created by 
the crystal lattice, which can result in the forbidden en- 
ergy bands, or gaps, for the electrons Q. Analo gous 
effects exist for electromagnetic waves in photonic crys- 
tal structures, where the refractive index is a periodic 
function of spatial coordinates, resulting in the photonic 
band gaps (PBGs) @,[|. 

Recently, spectacular progress has been achieved in 
cooling and trapping atoms in the optical lattices 
(OLs) — spatially periodic dipole potentials induced by 
off resonant laser fields [4i]. The relevant parameters of 
these systems can be controlled with very high precision 
and can be tuned to implement with unprecedented accu- 
racy some of the fundamental models of condensed mat- 
ter physics. A particularly relevant for the present stud- 
ies achievement has been the demonstration of the transi- 
tion from the superfluid to the Mott insulator phase with 
a commensurate number of bosonic atoms per lattice site 
i- 

As discussed below, OLs loaded with cold atoms in- 
teracting with a weak probe field and simultaneously 
driven by a strong coherent field can serve as a conve- 
nient tunable platform for the simulation and studies of 
light transmission and Bragg reflection in periodic media 
exhibiting PBGs. Earlier studies of related systems in- 
clude a weak Bragg reflection of the probe field from the 
sparsely occupied by atoms ID OL @, while more re- 
cently, optically induced ID PBGs in uniform (gaseous) 
atomic media were realized 0, Q employing the electro- 
magnetically induced transparency (EIT). EIT is a quan- 
tum interference effect characterized by the presence of a 
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frequency region with greatly reduced absorption accom- 
panied by steep dispersion for a weak probe field prop- 
agating in a three-level atomic medium whose adjacent 
transition is driven by a strong coherent field 0, [loj . 

Optical lattices loaded with cold atoms under the EIT 
conditions offer unique advantages over the previously 
studied schemes. First, the transparency bandwidth and 
the steep dispersion of the EIT resonance can be easily 
controlled by the corresponding driving field Sec- 
ond, for deep enough lattice potential, the Mott insulator 
regime can be reached [5] , in which each site will contain a 
single tightly localized atom. This will practically elimi- 
nate the coherence relaxation on the two-photon Raman 
transition caused by inter-atomic collisions and atomic 
time of flight, which result in a residual absorption of 
weak (quantum) fields and limit interaction times in most 
EIT experiments. Simultaneously, due to the tight lo- 
calization of the atoms, the local atomic density is in- 
creased by two to three orders of magnitude, as com- 
pared to the atomic density in the usual homogeneous 
(gaseous) media. This, in turn, strongly enhances the 
modulation amplitude of the effective refractive index, 
resulting in a strong Bragg reflection of the probe fleld. 
The high-contrast periodic index modulation is necessary 
for achieving 2D and 3D PBGs and strongly trapping the 
weak probe field. By contrast, in a uniform EIT medium 
7] , the dispersion modulation with an off-resonant stand- 
ing wave field can only yield the index contrast up to a 
few percent, which is sufficient for ID Bragg refiection 
over extended propagation lengths amounting to thou- 
sands of periods, but not enough for strongly trapping 
the light in 2D or 3D with several hundred periods in 
each direction. 

The paper is organized as follows. In Sec. |TT1 the po- 
larization of the atomic medium for a weak probe field 
under the EIT conditions is rigorously derived. This is 
then used to obtain the ID coupled-mode equations for 
the probe field. The solution of these equations and the 
resulting refiection, transmission and absorption spectra 
for the probe field in the cold atomic medium are pre- 
sented in Sec. IIIII A comparison with the case of a ther- 
mal atomic gas is made in Sec. IIVI The conclusions are 
summarized in Sec. IVl 
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FIG. 1: (a) Level scheme of cold atoms interacting with a 
weak probe field E and strong driving field with Rabi fre- 
quency Q.d- (b) The atoms are trapped in a tight-binding 
optical lattice potential Vi. 



II. MATHEMATICAL FORMALISM 



Consider the interaction of a weak probe field with an 
ensemble of cold bosonic (alkali) atoms with level con- 
figuration shown in Fig. [ija). The probe field E acts 
on the transition from the ground \g) to the excited |e) 
state. The transition |e) ^ k) is driven by a strong 
(near) resonant field with Rabi frequency ^d- The atoms 
in the lower meta-stable states \g) and \s) are subject 
to external spatially periodic potential Vg^s{^) schemat- 
ically shown in Fig. [IJb). The external potential Ve for 
the excited state atoms is assumed negligibly small. All 
the atoms are initially prepared in the ground state \g). 
At sufficiently low temperatures, only the lowest Bloch 
band of the lattice potential is populated. In the limit 
of large lattice depth, this band reduces to the lowest 
energy levels of individual potential wells. Then, for a 
lattice filling factor p ~ I, and provided that tunneling 
rate of the atoms between the lattice sites is smaller than 
on-site interaction energy due to the s-wave scattering, 
the Mott-insulator phase with single atom per lattice site 
will be attained 0, P 



A. Polarization of the Cold Atomic Medium 

The medium is described by three bosonic field op- 
erators ^/'i(r) representing atoms in the corresponding 
internal state \i) (i = g, e, s) and possessing the stan- 
dard bosonic commutation relations [tpi{r),^l,{r')] = 
S,,,S{r~r') and [Mr)Je{r')] = [i>l{r)Jl(r')] = 0. 
In general, the Hamiltonian of the system has the form 



H = Ha + H 



AA 



Hp + H 



AF- 



(1) 



Here Ha is the atomic Hamiltonian 



Ha = }_^J d'r^lir) -^V' + Viir) + hWi V^.(r), 

. (2) 

with M being the atomic mass and HcOi the internal en- 
ergy of the atoms in state The term Haa describes 
the interatomic collisions, which, under the conditions of 
single atom per lattice site and negligible intersite tun- 
nefing assumed above, does not play a role [1, Next, 
if the probe field is assumed quantized, its Hamiltonian 
Hp must be included in Eq. 11}. In the regime of a 
linear-response of the medium to a weak field studied 
here, the propagation dynamics of a quantized or a clas- 
sical field is the same and the probe field can be treated 
classically. Finally, Haf describes the interaction of the 
atoms with the probe E{r,t) = £(r)e~"^* and driving 
Ed{r,t) = f^e*^'''^-'-"''*' fields, 

Haf ^ j d\ijl{v)[-p,gE{v,t)\^g{v) 

d\ V-t (r) [ - p,sEd{r, t)] (r) + H. c, (3) 

where pu' is the dipole matrix element of the correspond- 
ing atomic transition \i) ^ and H. c. stands for 
the Hermite conjugate. With the Hamiltonian ([T]), the 
Heisenberg equations of motion for the atomic operators 
ipi read 

1 - 
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_j^£*(r)e-*4(r). 
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V's(r) = - 



i 

h 
-in* 



^s(r). 
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(4c) 



where fid = pes^d/^ is the drive Rabi frequency. 

The lattice potential is induced by far off-resonant 
standing-wave fields. It has the form Ki(r) = 
^Y^(^=x y z ^iii cos^(fcs^), where the modulation (ac Stark 
shift) amplitude iS'i;^ is proportional to the dynamic po- 
larizability of the corresponding atomic state \i) {i — 
g, s) and the intensity of the standing-wave fields with 
the wave vector kg. As stated above, in a deep OL po- 
tential, the atomic wavefunction is tightly localized at 
each lattice site j. The atomic field operators i/'i(r) can 
then be expanded in terms of the real, normalized (Wan- 

nier) functions ^^^"'•'(r) = ^^(r — r^), which in the vicinity 
of r ~ Yj satisfy 



fi2 

2M ^ > 



{v)^hu,w'f\v), (5) 
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where hPi is the energy of the lowest vibrational state of 
the nearly harmonic potential well, while w'f^ (r) can well 
be approximated by 3D Gaussian functions 

3/4 



«;p)(r) 



1 



exp 



(r- 



^2^ 



(6) 



centered around with the width (5?'. For tight lo- 
calization Sr <C Tr/ks, the Wannier functions pertain- 
ing to different lattice sites have negligible overlap, 
/ cfir w^/^ (r) tup (r) = Sjj> . The field operators for 
the lower atomic states can then be decomposed as 
M^^t) = Ej4'\r)g,(t)e-*('^«+'^«)* and ^s{r,t) - 

X;jwi-'\r)sj(t)e-*('^s )*e"*('^~'^'')*, where gj{t) and 
Sj{t) are the slowly- varying in time annihilation oper- 
ators for the bosonic atoms at site j. Since the atoms 
in the excited state |e) are assumed free, Ve — 0, the 
corresponding atomic field operator can be expanded 
as 4(r,i) = Ek,^k,(r)ek,(Oe-*('^«+'^«)*e-^'"*, where 
Mk, (r) — (L)^'^/^e'''' "" are the plane waves within the 
quantization (medium) volume i"^, which satisfy 



d3r<,(r) Uk;(r) = 5; 



kik; 



mi 
2M ' 



(7) 
(8) 



and Ck, {t) are the slowly- varying in time annihilation op- 
erators for the corresponding mode. 

Substituting the above decompositions of the atomic 
field operators into Eqs. ^ yields 
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3 



'^^*(i")5I"k,(r)ek, +^3, (9a) 



^ Mk, (r) ek, 

k. 



= ^ ^^^^ ^ (r) Ck, 

k, 

-f*^£(r)^z««(r)5^. 

3 

e*''-'- wi^'> (r) S, + A, (9b) 



Ywy\r)§, = (»A^-7.)^^?')(r)s, 

3 3 

+ifiSe-^''''-"-^Uk,(r) Ck, +.F„(9c) 

ki 

where A = uj + ajg + iig — LOe and An = uj + LOg + Vg— tus- 
h's — are, respectively, the one-photon and two-photon 
(Raman) detunings of the probe field [cf. Fig.[IJa)], while 
7i are the phenomcnological half-decay rates of the cor- 
responding atomic states \i) [i — g,e,s). Note that 
7e 10'' s~^ is large, determined by the short radia- 
tive lifetime of the excited atomic state, while for iso- 
lated cold atoms trapped in a conservative OL poten- 
tial, 7g.s ^ 1 are negligibly small, due to the absence 



of interatomic collisions and atomic free-flight or diffu- 
sion away for the interaction region. Finally, J^i are the 
(5-correlated noise operators associated with the relax- 
ation. These noise operators, required to preserve the 
formal consistency of the equations for the atomic oper- 
ators, will not play a role in the following semiclassical 
treatment of the field propagation and will be dropped 
from now on. 

Under the weak probe approximation, the initial 
atomic population in the ground state \g) remains practi- 
cally undepleted and Eqs. ([9]) can be solved to the lowest 
(first) order in £{r). In the stationary regime, Eq. (Pc]) 
yields 



3 



(r) 



-ikd-r 



Ek, "k, (r) Ck, 



iA^ - 7^ 



Substituting this into Eq. Ij9bp . multiplying all terms by 
J d?r u^^ (r) and using Eq. ([5]) leads to 



ek, 



.Peg EJd^rulir)£ir)w^/\r)g, 



7e + 



iAR-7s 



The integral in the above equation can be evaluated as- 
suming that the localization width 6r of w'^J^ (r) is small 
compared to the wavelength A = ^-ncjio of the probe field. 
Therefore f (r) changes little in the vicinity of r ~ Vj and 
can be taken constant in this region. The remaining in- 
tegral is then given by 



7-0) _ 
^k, - 



/ 1 \ 



The polarization P{v,t) — 7'(r)e^*'^* of the atomic 
medium, induced by the applied probe field E(r,t), 
is given by the expectation value of the atomic 
dipole moment at position r and time t, P(r, t) = 
(■01 (r, t) pge V'e(i", t))- With the above decompositions of 
the atomic field operators i^g^e and using Eq. PH)) . the 
medium polarization V(r) reads 



Viv) = p,,(5]i«(^-)(r)5]5]«k,(r)ek, 



■J2w(g^\r)£{ry){g]g,,) 



E 



"k,(r)X;^ 



k, ^A-7e-r^^_^^ 



(12) 



Here the sum over the plane-wave modes k; should be re- 
placed by an integral according to ~^ i^/'^'^)^ I d^^i 
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m 



which, upon substituting •* from Eq. pTjl . reads 



LV f 1 



2tt J \TTdr^ 



3/4 



2tt Sr 



X / (fki 



iA - 7e 



(13) 



It is easy to see that if the atomic locahzation width is not 

too small, 6r^ > 6rl,^ = h/M-fe 0, J^fe, = ^ in the 
denominator of the above equation contributes little to 
the integral and can therefore be dropped. The integral 
is then easily evaluated, and Eq. reduces to 



field can be expanded in terms of the free-space nor- 
mal modes as E{z,t) = Y.k^k{z)e'^^''-'^*\ where £k{z) 
are now slowly varying in space amplitudes of the corre- 
sponding field modes k = -iiLo/c. The medium polariza- 
tion can be expressed as P[z, t) = eox('^; z)E(z, t), where 
x(i^; z) is the linear susceptibility, whose z-dependence is 
due to the atomic density g{z). This ID density is ob- 
tained by averaging g(r) over the transverse to z direc- 
tions of the cubic lattice with primitive vector A = Tr/fc^, 



I M/a /-A/a 

J-A/2 J-A/2 



dyg{r). 



(16) 



,(/) 



(r) 



Since g(z) = g{z + A) is a periodic function, it can be 
expanded in a Fourier series 



iA - 7e 



e(z) = ^fte^'^^ g 



2tt 

T 



(17) 



Since the Wannier functions w)j" {r) and Wg'' \r) pertain- 
ing to different lattice sites j ^ f have negligible overlap, 
the final expression for the medium polarization takes the 
form 



with the expansion coefficients 
1 /-A/z 



Qi 



A 



A/2 



dz g{z) e-'^s" 



V(r) 



I peg I 



g{r)£{r) 



iA- 



(14) 



where g{r) = J^j \'^g'\^)\'^ (lijdj) is the spatially- 
periodic density of the medium. According to the as- 
sumption above, the occupation number of each lattice 
site is (gjgj) = 1. Clearly, the spectral response of the 
medium coincides with that of the conventional EIT with 
homogeneous atomic ensemble 0, [l3|, but is spatially 
modulated by the periodic atomic density, which leads 
to profound consequences for the probe transmission and 
reflection, as discussed below. 



^ MM 

A3 

where, according to Eq. 

|u;,(z)P = 



A/2 



A/2 



1 



dz\wg{z)\^e''^s'', (18) 



■ exp 



z 



(19) 



is a normalized Gaussian function. Thus, the zeroth 
Fourier component go — {g^jgj)/A^ (Sr <C A) is the 
spatially averaged (bulk) atomic density. Higher order 
Fourier components can be expressed as gi = qqUi, where 



A/2 



A/2 



dz\wg{z)\^e-'^s''. 



(20) 



B. Coupled-Mode Equations for the Probe Field tt • ti^ mn. ii j- i • j^- j 

^ Usmg Eq. p4|) . the medium polarization then reads 



The propagation of electromagnetic field in the 
medium is governed by the Maxwell wave equation 



P{z, t) = eo x(w) £k' {z) e 

I k' 

where 

X(w) = XEIT(t^) 



2 1 



d_ 



E{r,t) 



92 

Mo^P(M). 



(15) 



i{lg~\-k )z —iut 



In free space, P = 0, its general solution for a monochro- 
matic field of frequency uj has the form E{r,t) — 
'^ke*^'' ''~'^*\ where £k is the amplitude of the field 
mode with wave vector k satisfying |k| = w/c, with 
c = (/xqEo)"^/^ being the speed of light in vacuum. 

Consider a stationary propagation of the probe field in 
the periodic atomic medium whose spatial and spectral 
properties are characterized by Eq. (|14p . In particular, 
assume that a circularly (cr^) polarized probe field propa- 
gates along the z axis (taken to be the quantization axis), 
while a linearly (tt) polarized driving field, propagating 
along the x axis (k^-Lz ), uniformly irradiates the whole 
atomic sample. Following the approach of [3|, the probe 
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UJ 



«7e 



7e 



iA- 



(21) 



(22) 



is the usual EIT susceptibility of an atomic medium with 
uniform density go and resonant absorption cross-section 



I p. 



"eg 



/{2eoch'-fe) By definition, the res- 

onant amplitude absorption coefficient in the two-level 
atomic medium {fid = 0) is oq = aogo- 

Using Eq. (I2ip . Maxwell's equation yields 



eg I 



2i^k 



Akz 



I k' 
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which, upon expanding the sums over k, k' = ±uj/c, leads 
to the following coupled equations for the amplitudes 
£± = S±uj/c of the forward (+) and backward (— ) prop- 
agating modes, 



—£+ = ia(w) ^ [£+ e*'9^ + £_ e 



i{lg — 2uj I c)z'\ 



(24a) 



dz 



where 



a(t^) = ^ 

c 2 



7e - i A + 



7s— iA_n. 



(25) 



Note that k/ decrease as |/| increase. Thus, apart from 
Ko — 1, K±i are the largest coefficients given by {6r <C A) 



1 



dze"^'/'^'''cos(27rz/A). (26) 



Assuming, therefore, that lo/c^ k^ and keeping only the 
terms satisfying the longitudinal phase matching condi- 
tion (i.e., slowly oscillating in space), Eqs. ([24|) reduce to 



—£+ = ^a(^J)f+ +^a(a;)A^lf_e-2»('^/'=-'=»)^ (27a) 
dz 

^ -ia{uj)£--ia{Lj)K-i£+e^''-'^/^-'''^\ (27b) 

In these equations, the first term describes the usual EIT 
absorption and dispersion of the probe field propagating 
in the atomic medium, while the second term is respon- 
sible for the coupling of the forward and backward prop- 
agating modes mediated by the medium periodicity. 



III. REFLECTION, TRANSMISSION AND 
ABSORPTION OF THE PROBE 



The coupled mode equations ([27)1 fully determine the 
optical properties of the system under consideration. The 
solution of Eqs. ([?7|) for the boundary value problem de- 
fined through £+{0) = £in and £-{L) = is given by 



£+{z) = e-*('^/'=-'==)^ 

scosh [s{L — z)] — i5f3 sinh [s{L — z)] 
s cosh [sL] — iS/3 sinh [sL] 

ia{uj) K_i sinh \s{L — z)] 
s cosh [sL] — iS/3 sinh [sL] ' 

where the coefficients 

a(uj) H ks, 



,(28a) 



(28b) 
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FIG. 2: Absorption (blue dotted lines) and dispersion (black 
solid lines) spectra for the probe field in a uniform EIT atomic 
medium. The Rabi frequency of the resonant driving field is 
(a) Qd = 27e, and (b) = 57e. 



quantify, respectively, the phase mismatch and the cou- 
pling between the forward and backward propagating 
modes. 

According to the Bloch theorem 0, Q , the general so- 
lution for the electromagnetic field propagation in the pe- 



riodic medium can be cast as £{z) 



iKi 



£k{z), where 



if is a propagation constant known as the Bloch wave 
vector, while £k{z) — £k{z + A) is a spatially periodic 
function. In a purely dispersive (non-absorbing) medium, 
the real and imaginary parts of K describe, respectively, 
the spatial phase and the attenuation (due to Bragg re- 
flection) of the field upon propagation. It follows from 
Eqs. (|28p that the corresponding propagation constant 
for the monochromatic probe field of frequency uj is given 
by 



K 



IS 



kg + i\l [c({uj)ki]^ — a{uj) -\ kg 



(29) 



In the present situation, the periodic atomic medium ex- 
hibits both strong dispersion and absorption character- 
ized by the EIT polarizability a{uj) (see Fig. [2|). There- 
fore, in the dispersion relation of Eq. (P^ . both Bragg re- 
flection and medium absorption would contribute to the 
imaginary part of K. A more quantitative and experi- 
mentally accessible characterization of the optical prop- 
erties of the system is provided by the reflection R, trans- 
mission T and absorption A coefficients for the probe 
field, which are defined through 



A{u;) 



£_(0) 



£+{0) 
£+{L) 



1 



£+{0) 



(30a) 

(30b) 

(30c) 



The absorption and dispersion spectra for the probe 
field in a uniform atomic medium subject to a resonant 
drive LOd = (^es — i^s is shown in Fig. [2l The strong driv- 
ing field with Rabi frequency fl^ > 7e splits the familiar 
Lorentzian absorption line into the Autler-Townes dou- 
blet with the two peaks separated by 2ild, while at the 
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line center the medium becomes transparent to the probe 
field. This EIT is accompanied by a steep normal dis- 
persion, resulting in the greatly reduced group velocity 



M 



1 



aole 



< c. 



(31) 



The principal aim here is to explore the possibility of 
a PBG for the probe field in the frequency region of the 
reduced absorption. Assuming resonant driving field, so 
that A/j = A, in the vicinity of EIT resonance |A| < fid, 
the polarizability of Eq. ([25|) reduces to a(w) ~ A/wg 
and the dispersion relation (P^)) can be approximated by 



"A ■ 

— Kl 

.^9 . 


2 


'A 




2 






c 





(32) 



where Sg = kgC — ojeg- The Bloch wave vector K has 
an imaginary part when the term under the square-root 
of this equation is positive. The corresponding range of 
frequencies of the probe field is 



(1 ± Ki)c 



< A < 



(lT'«i)c' 



for 5s ^ 0. 



(33) 



Thus, a PBG of width ^Wgap = \^s\^y::^, 



centered at 



A 



ap — i-i^i ' formed on the positive (for 5s > 0) 
or negative (for 5s < 0) side of the EIT resonance. The 
peak value of the imaginary part of which determines 
the maximal reflectivity at Agap, is given by 



max(Im K) 



\5s 



(34) 



These conclusions are justified upon requiring that 
the PBG lies within the EIT window Scotw = 
|^^dP/(7e\/2aoL) [13, where the absorption is small. 
This leads to the following condition on the lattice 
wavevector fc,, 



Ifc. 



_ \6s 



< 



(1 




(35) 



with which max(Imiir) < Kiy/ao{l — Ki)/2L. With 
the system parameters listed in [ll|, and for an atomic 
medium of length L ~ 200 /xm (500 lattice periods A), 
or optical depth 2aaL — 100, the above condition is sat- 
isfied for \ks — ijJeg/c\ ~ 4.5 X lO'^ rad/m or |(5s|/27r ~ 
2.16 X 10^^ s~^. Physically, such a large difference be- 
tween the lattice wave vector ks and that of the probe 
field k ~ LOeg I c stems from the need to satisfy the phase 
matching condition J/3 ~ [while a{Lo)K\ ^ 0] in strongly 
dispersive EIT medium. 

The complete dispersion relation of Eq. ((29|). with 
Kl ~ 0.9 obtained for 5r ~ A/10 [ll|, is plotted in 
Fig. [3] for the two values of il^ used in Fig. [51 As noted 
above, the imaginary part of K describes exponential 
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FIG. 3: Imaginary (blue dotted lines) and real (black solid 
lined) parts of iC' as a function of probe detuning A (resonant 
drive, A = Aa), for ki = 0.9 (5r ~ A/10), 5s = 7.2 x 10* 7e, 
and (a) Q.d = 2^e and (b) Q.d ~ 5^e- Insets magnify the 
important frequency regions within the EIT window. 



attenuation of the forward propagating probe field due 
to the back reflection and medium absorption. Within 
the transparency window, |A| < <5wtw, the absorption is 
much smaller than the spatially-periodic dispersion, and 
the bulges of ImK in the insets of Fig. [3] signify the ap- 
pearance of the PBG. The band edges given by Eq. ((33l) 
can be tuned by changing the drive field intensity, since 

Vg cx irirfp. 

Figure |4] shows the refiection, transmission and absorp- 
tion spectra for the probe field for the two values of Sid 
used in Figs. [2] and [3l The peak refiectivity for the PBG 
within the EIT window is about 90%, limited mainly by 
absorption. Choosing smaller values of 5s, thereby mov- 
ing the PBG closer to the EIT line center, and taking si- 
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FIG. 4: Reflection R, transmission T and absorption A spec- 
tra for the probe field in a cold EIT atomic medium loaded 
into an OL. The medium length L = 500 A ~ 200 /xm, or 
optical depth 2aoL = 100, and all the other parameters are 
same as in Fig. |3ja) and (b) , respectively. Insets magnify the 
important frequency regions within the EIT window. 
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FIG. 5: Diagram illustrating the appearance of the lower Sujl 
and upper 5u>u PBGs far away from the EIT resonance, where 
the atoms behave as off-resonant two-level systems with neg- 
ligible absorption. 



multaneously longer medium, to compensate for smaller 
values of max(Im K) as per Eq. (|34|) , will yield even larger 
values of the maximum reflection (up to 98%), at the ex- 
pense of narrower PBG. 

In the spectral region far away from the atomic line 
center and EIT features, there is another broad PBG 
seen in the reflection spectrum of Fig.[4l Its existence was 
predicted and studied in for the limiting case of ki — > 
1, corresponding to the atomic localization width Sr — s- 
0. Briefly, in the frequency regions |A| ^ je,^d, the 
polarizability of Eq. (j25[) can be approximated as a{Lo) ~ 
— ao7e/A, which describes the response of off-resonant 
two-level atoms, and the dispersion relation reads 




A 



S..-A 



(36) 



The Bloch wave vector acquires imaginary part when 
ao7eC(l — Ki) < A(A — 6s) < ao7ec(l -I- ki), which is 
graphically illustrated in Fig. [3 Two PBGs appear, one 
on the red A < and the other on the blue A > sides 
of the atomic resonance. The corresponding ranges of 
frequencies for the lower (L) and upper (U) gaps are 



-D+< Al < hSs - 



D^<Au< Us 



D 



+ 1 



(37a) 
(37b) 



where D± = ^ {Ss/'^Y + ao7ec(l ± ki), and the widths 
of the gaps are 5ujl,u = D+ — When (5s /2)^ > 

ao7eC, which for 5s chosen as in Eq. (|35p is certainly sat- 
isfied, D± can be expanded as D± ~ |5s|/2 -f ao7ec(l ± 
Ki)/|(5s|. For 5s > 0, the lower gap — ao7ec(l-l-Ki)/|5s| < 
Al < — ao7ec(l — Ki)/|(5s| is closer to the atomic reso- 
nance, and the corresponding reflection is seen in Fig. U 
the upper gap is beyond the range of frequencies spanned 
in that figure. Note finally that the reflection, transmis- 
sion and absorption spectra for 5s > and 5s < are 
related by mirror symmetry about A = 0. 



IV. THERMAL ATOMIC MEDIUM 

For the sake of comparison, it is instructive to con- 
sider the case of a thermal atomic gas under the other- 
wise similar conditions. When the thermal energy k^T is 
much larger than the depth of the OL potential hSi , the 
atoms move freely and the atomic density g is uniform 
[T^ . However, the off-resonant standing wave field with 
wave-vector ks results in a spatially-periodic ac Stark- 
shift hSi cos'^{ksz) of the lower atomic levels \i) {i = g, s). 



The resulting susceptibility 



^ 7e - iA{z) + 



(38) 



7s-iA_R(z) 

is a periodic function of z, since the corresponding one- 
and two-photon detunings A(z) = A'+^Sg cos{2ksz) and 
Ar{z) — A^ + ^Sgs cos{2ksz) are ac Stark modulated. 



Here A' = cj — < 



iSg and A'^ 



UJ-UJd + UJns + 



are the mean detunings, while Sgs = Sg — Ss is the dif- 
ference of the ac Stark modulation amplitudes for levels 
\g) and \s). The decay rate 7s describes the Raman co- 
herence relaxation, which is now affected by the atomic 
motion and collisions; its typical value is in the range of 
7s ~ 10'^ — 10''s~^ 0], which is still much smaller than 7e 
but not negligible. The effect of the spatially-periodic 
modulation of the Raman detuning Afl(z) is to peri- 
odically shift the EIT spectrum for the probe field (see 
Fig. m) , which under appropriate conditions discussed be- 
low can result in a PBG, as was shown in Q. 

Substituting the susceptibility of Eq. (|38p into the 
Maxwell equation (fT5)l with P{z,t) — eoxi^^', z)E{z,t) 
and E{z,t) = J2k^k{z)e'^''''-'^^\ under the EIT condi- 
tions l^ldP ^ lejs, {jSg,s)'^ and the longitudinal phase 
matching kg ~ w/c, the following coupled mode equa- 
tions for the forward and backward propagating modes 
k — ztoj/c of the probe field are obtained, 

= ta'{uj)£++ir,{uj)£^e~^'^^^^''''^', (39a) 

d 
dz 
where 



-ia'{uj)£_ - iT]{(^)£+ e^^C'^/^-'^-)^, (39b) 



a'(w) = 



ao7e 



:{A'^i\n,f-A'A'^) 
+i[\nd\H+leiA'^ 



1 o2 
S'-'gs 



)]},(40a 



ri{uj) 



ao7e 



A'Ay%7?A'J 
x^{\na\^ Sgs -A'j,{2A' Sgs 

+2ijeA'j^Sgs)y 
The solution of Eqs. ([55]) is given by Eqs. 
a'{uj) + ''' 



A'k5,) 



(40b) 



with 



5p 



ks ) 
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FIG. 6: Imaginary (blue dotted lines) and real (black solid 
lines) parts of /C" as a function of probe detuning A' (reso- 
nant drive, A' — A^), for fla — 2'ye, 7s = 10~*7e and Ss = 0. 
The ac Stark shift amplitudes are (a) jSgs = 0.27e, and 
(b) jSgs ~ 0.47e. Insets magnify the important frequency 
regions within the EIT window. 
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The corresponding dispersion relation K = + is, with 
a resonant driving field ujd = oj^g -~ (A'^ — A') and 
kg ~ Weg/c, is shown in Fig. [51 where nonzero values of 
liaK within the EIT window signify the appearance of 
a PBG for the probe field. 

In the vicinity of EIT resonance |A'| < ild (assuming 
resonant drive, A' — A^), Eqs. (fiO]) reduce to a'{uj) ~ 
A'/wg and r]{Lu) ~ Sgs/{^Vg), and the dispersion relation 
can be approximated as 



Sgs 


2 


'A' s; 


.4W3. 







(41) 



For = 0, a PBG for the probe field with frequencies in 
the range |A'| < \ \Sgs\ is formed ((Jwgap = \\Sgs\)- The 
peak value of Im K attained at the gap center Ag^^p = 
is max(Imi4:) = \Sgs\/ {'ivg). Note that a PBG for the 
probe field exists only when Sgs = Sg — Ss 0, i.e., the 
ac Stark shifts Sg and Ss of the lower atomic levels \g) 
and \s) induced by an off- resonant standing wave field 
should be different. It is clear that in order to minimize 
absorption, the ac Stark modulation of the Raman res- 
onance should be accommodated with the EIT window, 
^I'S'gsl ^ ^"^tw 0- This, in turn, restricts the peak reflec- 
tivity of the medium to max(Im K) < ^Joq/WL. 

Figure [7] shows the reflection, transmission and absorp- 
tion spectra for the probe field for the two values of Sgs 
used in Fig.[Sl In Fig.[7]Ja) with \Sgs = 0.27e, the peak of 
the reflection coefficient i? at A' = is about 80% and the 
absorption A is 18%. Even though increasing Sgs results 
in broader band gap and larger values of IthK, as seen 
in Fig. [S]Jb) with \Sgs = 0.47e, the corresponding reflec- 
tion coefficient of Fig. [Tl^b) is smaller, R ~ 67%, due to 
the increased absorption A ~ 33%. Physically, this can 
be understood by recalling (see Fig. [2] and [3, [ifl]) that 
in the vicinity of EIT resonance, the dispersion (which 
determines IuyK) scales linearly with detuning A^ (i.e., 
with Sgs), while absorption scales quadratically with Sgs- 
Therefore, by choosing smaller values of the ac Stark 
modulation amplitudes Sgs, and taking simultaneously 



FIG. 7: Reflection R, transmission T and absorption A spec- 
tra for the probe field in a thermal EIT atomic medium sub- 
ject to spatially periodic ac Stark modulation of the Raman 
resonance. The medium length is L ~ 200 /im (optical depth 
2aoL = 100), and all the other parameters are same as in 
Fig. |6ja) and (b), respectively. Insets magnify the important 
frequency regions within the EIT window. 



longer medium, to compensate for reduced lu\K, will 
yield larger values of the reflection coefficient at the cen- 
ter of the PBG which will however shrink in width. 



CONCLUSIONS AND OUTLOOK 



Equation or (PS)) characterize the spectral prop- 
erties of clcctromagnetically induced transparency (EIT) 
in an atomic medium, as illustrated in Fig. O Looking 
at that figure, one could come up with three possible 
ways for achieving spatially-periodic index modulation 
5n{uj; z) ~ iRex(w; z) resulting in a photonic band gap 
(PBG) in the vicinity of EIT resonance |A^| < r2rf,7e, 
where the susceptibility is approximately given by 



x(w;2:) 



C crog7e . 
^ ~ lO 12 



(i) To periodically modulate the medium density g = q{z) 
by trapping cold atoms in an optical lattice, as proposed 
and studied here, (ii) To periodically shift the EIT spec- 
trum A/j — /S.fi{z) in a thermal atomic ensemble by us- 
ing spatially periodic ac Stark shift of the Raman transi- 
tion induced by off-resonant standing-wave field, as pro- 
posed in [7] and reviewed in Sec. IIVI above, (iii) To use a 
standing-wave drive field U,d = ^d{z) resulting in a spa- 
tially periodic modulation of the bandwidth of the EIT 
and the associated with it dispersion slope, as studied in 
i. 

The novel scheme proposed here possesses a number 
of advantageous properties. These include: (a) Greatly 
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reduced coherence relaxation on the two-photon Raman 
transition due to the eUmination of atomic diffusion and 
colhsions causing probe absorption, (b) Simultaneous en- 
hancement of the refractive index modulation and the 
resulting Bragg reflection due to the tight localization 
of the atoms, (c) Tunability of the position and the 
width of the PBG within and beyond the EIT resonance. 
These properties can be employed for achieving very effi- 
cient nonlinear interactions between weak quantum fields 
and realization of deterministic quantum logic with single 
photons, as described in |14|. 

Finally, it would be interesting to explore the possi- 
bility of achieving a 2D PBGs for the probe field in the 
medium of cold atoms trapped in an optical lattice. The 
ultimate goal would be to explore structures with appro- 
priately engineered defects which may allow for strong 
localization and waveguiding of light. Such defects can 
easily be implemented in the laboratory experiments by 
simply focusing a resonant laser onto the desired lattice 
sites, which will release (evaporate) the corresponding 



atoms from the trap. In contrast, the microfabrication of 
desired defects in solid-state photonic crystal structures 
involves complicated and time consuming growth and 
lithographic techniques which are often not easily recon- 
figurable and reproducible. Thus tunable photonic band 
gaps in optical lattices may serve to experimentally test 
various proposed structures in a simple, quickly repro- 
ducible way. 
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